'''
例如传染病预测模型，经济增长预测模型，兰彻斯特战争预测模型，
药物在体内分布与排除预测模型，人口预测模型，烟雾扩散与消失预测模型

例题：美日硫磺岛战役模型
A(t) J(t)表示美军 日军第t天的人数，u(t)表示美军增援
    dA(t)/dt=-aJ(t)+u(t)
    dJ(t)/dt=-bA(t)
    A(0)=0 J(0)=21500 J(3)=0

    两边积分以求和来近似代替，得
    A(t)-A(0)=a*sum(J(t))+sum(u(t))      1式
    J(t)-J(0)=-b*sum(A(t))               2时
    由实际数据得到sum(A(t))=20370000，代入得到b=0.0106
    b代回2式可以得到J(t) t=1,2,...,36
    然后令t=36 代入1式估计a，得到a=0.0544
    a代回1式可以得到A(t) t=1,2,...,36

    u(t)=54000 0<=t<1
         6000  2<=t<3
         13000 5<=t<6
         0     others
'''
a=0.0544
b=0.0106
sumA=2037000
sumJ=372500
A=[]
A[0]=0
def uFun(t):
    if 0<=t<1:
        return 54000
    elif 2<=t<3:
        return 6500
    elif 5<=t<6:
        return 13000
    else:
        return 0
for i in range(1,37):
    A.append(A[0]-a*sumJ+uFun(i))
J=[]
J[0]=21500
# 数据不足，代码省略
